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ABSTRACT 

We present high S/N, moderate resolution (R « 2000) near- infrared spectra, as well 
as 10 fim imaging, for the brightest members of the central stellar cluster in the W40 
H II region, obtained using the SpeX and MIRSI instruments at NASA's Infrared Tele- 
scope Facility (IRTF). Using these observations combined with archival Spitzer Space 
Telescope data, we have determined the spectral classifications, extinction, distances, 
and spectral energy distributions for the brightest members of the cluster. Of the eight 
objects observed, we identify four main sequence (MS) OB stars (one late-O, 3 early-B), 
two Herbig Ae/Be stars, and two low-mass young stellar objects (Class II). Strong He I 
absorption at 1.083 /im in the MS star spectra strongly suggests that at least some of 
these sources are in fact close binaries. Two out of the four MS stars also show sig- 
nificant infrared excesses typical of circumstellar disks. Extinctions and distances were 
determined for each MS star by fitting model stellar atmospheres to the SEDs. We 
estimate a distance to the cluster of between 455 and 535 pc, which agrees well with 
earlier (but far less precise) distance estimates. We conclude that the late-0 star we 
identify is the dominant source of LyC luminosity needed to power the W40 H II region 
and is the likely source of the stellar wind that has blown a large (~ 4 pc) pinched-waist 
bubble observed in wide field mid-IR images. We also suggest that 3.6 cm radio emis- 
sion observed from some of the sources in the cluster is likely not due to emission from 
ultra-compact H II regions, as suggested in other work, due to size constraints based on 
our derived distance to the cluster. Finally, we also present a discussion of the curious 
source IRS 3A, which has a very strong mid-IR excess (despite its B3 MS classification) 
and appears to be embedded in a dusty envelope roughly 2700 AU in size. 

Subject headings: Infrared: stars — ISM: individual objects (W40[W 40]) — circum- 
stellar matter — Stars: early-type — Stars: pre-main sequence 



Introduction 



W40 is a thermal radio continuum source associated with the molecular cloud G 28.74+3.52 in 
the Aquila Rift QWesterhout||1958[ |Zeilik fc Lada||1978[ ). A faint, patchy optical H II region (Sh 64) 



is associated with the radio emission (Sharpless 1959). Radio observations imply a total Lyman 



continuum luminosity of ~ 1.5 x 10 48 photons s _1 ( Goss k, Shaver|[l970 Smith et"aL"|l985 ) — about 
10% of the ionizing flux found in the Trapezium Cluster in Orion — corresponding to several early 
B stars or a single 09V. An IRAS source, cold ammonia core, and a number of mm-wave sources 



are also associated with the cluster (Molinari et al. 1996; Maury et al. 2011), suggesting on-going 



star formation from dense molecular material. In addition, Crutcher & Chu ( 1982 ) and Wu et al 



(1996) have found multiple CO line components and/or extended line wings toward W40 which 
suggests outflows driven by young stellar objects (YSOs). 
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Early infrared maps of the central 3' of W40 revealed 7 bright sources with optical counter- 
parts behind 9 — 10 magnitudes of visual extinction (Crutcher & Chu 1982 Smith et al. 1985) 



Spectral energy distributions (SEDs) from IR through millimeter wavelengths suggest substantial 



circumstellar material around most of these bright sources ( Smith et al. 1985 ; Vallee & MacLeod 



1994). 2MASS images of the region reveal a cluster of near-IR sources within the central 5' of the 



brightest sources (Fig. [I]). Recent radio observations at 3.6 cm show a cluster of compact sources 
in the central portion of the W40 IR cluster, many of which correspond to the known infrared 



sources (Rodriguez et al. 2010). X-Ray observations of the W40 cluster with the Chandra X-Ray 



Observatory (CXO) reveal approximately 200 sources associated with the cluster, the majority of 



which are thought to be low- mass young stellar objects (Kuhn et al. 2010). TheX-Ray luminosity 



function (XLF) obtained by the Chandra Orion Ultra-deep Project (COUP) indicates a complete- 
ness limit of 70% (down to 0.1 M ) and suggests a total cluster population of «600 objects (Kuhn 



etH|2010| >. UKIRT near-IR observations presented in the same study show that a significant frac- 
tion of the X-ray sources also have K-band excesses indicative of circumstellar disks and suggest an 
age of < 1 Myr ( |Kuhn et al.|2010[ ) . Herschel SPIRE and PACS far-IR maps of the entire Aquila Rift 
region (which includes the Serpens star forming regions and MWC297/Sh-62) reveal ~201 sources, 



most of which appear to be associated with W40 (Bontemps et al.|2010). Analysis of the far-IR 



colors suggest that 45-60 of the detected sources are deeply embedded Class objects (Bontemps 



et al. 2010) 



Depending on the technique used, the distance to W40 ranges from 300 to 900 pc (Rodney 



& Reipurth 2008, and references therein), making W40 one of only a few H II regions within a 
kiloparsec of the Sun. XLF fitting by iKuhn et al. (2010) yields a best-fit distance of « 600 pc, 



while Bontemps et al. (2010) adopt a much closer distance of 260 pc based on colocation of W40 



with the Serpens region on the sky. The distance to W40 will be discussed further in Section 6.1 



Many questions remain regarding the initial mass function (IMF), age, and star formation 
history of W40. Even the brightest members of the cluster have not been studied in any detail. In 
this paper we present infrared observations of these bright members of the W40 IR cluster in an 
effort to better constrain some of their basic properties, including spectral type, age, and distance. 
This is just a first step in better understanding the cluster and star forming region as a whole. 



2. Source Identification 



In Figure [T] we present the 2MASS K-band image of the central 6' of the cluster with the 



brightest sources identified and labelled using the Smith et al. ( 1985 ) naming convention. Inspection 



of 2MASS K-band image reveals that some of the original IR sources detected by Smith et al. ( 1985) 



are in fact multiples. In addition, there is another bright K-band source approximately 2' to the 
west of IRS 3A that corresponds to a bright star on the DSS R plate but is not reported by either 



Zeilik & Lada ( 1978 ) or Smith et al. ( 1985 ). In keeping with the current naming convention, we will 



refer to this source as IR Source 5 (Optical source 4A was not detected in the IR by Smith et al 
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(1985)). IRS 5 is surrounded by patchy nebulosity on both the DSS R plate and in the 2MASS K 



images. 



3. Observations 

3.1. SpeX Observations 

Observations of the central portion of the W40 IR cluster were carried out using NASA's 
Infrared Telescope Facility (IRTF) with the SpeX instrument in 2006 and then again in 2010. SpeX 
is the facility near-IR spectrograph for the IRTF and provides spectral resolving power (R) of 1000 - 
2500 across 0.8-2.4, 2.0-4.1, and 2.3-5.5 /im, using prism cross-dispersers with a Raytheon Aladdin 



3 1024x1024 InSb array (Rayner et al. 2003). SpeX also contains an infrared slit- viewer /guider 



covering a 60x60" field-of-view with a Raytheon Aladdin 2 512x512 InSb array. 

We obtained both images and spectra of the central region of W40 on 8 July 2006 under clear 
skies with K-band seeing of / /4-0'. / 6 (FWHM). Images of the central region of W40 were taken 
using the SpeX guide camera with a 9-point dither pattern centered on IRS 1A and 2A respectively. 
Dithered image sets were reduced using standard techniques and then mosaiced together; the final 
mosiaced image is shown in Figure [2] IRS 1 A is clearly resolved into two bright sources ( "North" 
and "South"), separated by « 2" along a roughly N-S position angle, and 3 additional weaker 
sources to the NE and NW. IRS ID is resolved into a small cluster of at least 7 distinct sources. 
IRS 2A appears to have a faint companion to the East as well (see Figure inset), but higher 
resolution observations are required for confirmation. 

Moderate resolution spectra (R~2000) were obtained for each of the brightest near-IR sources 
(Table [I]). The observations were acquired in "pair mode", in which the object was observed at two 
separate positions ("A" and "B") along the 15"-long slit. The slit width was set to 0'. / 3, which yields 
a nominal resolving power of 2000 for the short-wavelength cross-dispersed (SXD) spectra and 2500 
for the long- wavelength cross-dispersed (LXD2.1) spectra. The slit was set to the parallactic angle 
during the observations. In SXD mode, we used exposure times of 60 s for the brightest sources 
and 120 s for the fainter sources and cycled through the AB pairs twice yielding total on-source 
exposure times of 240 and 480 s respectively. This mode yields spectra spanning the wavelength 
range 0.8 -2.5 fim divided into 6 spectral orders. In LXD mode, we observed only the brightest 
targets with 30 s integration times, cycling through the AB pairs twice for IRS IB and 2A, and 
6 times for IRS 3A. This mode yields spectra covering 2.2—5.1 in 6 spectral orders. The 
airmass was 1.1-1.6 for all SXD observations and 2.0-2.3 for the LXD observations. The AO V 
stars HD 171149 and HD 169961 were observed throughout the night as telluric and flux calibration 
standards. The airmass difference between the objects and the standards was always < 0.2. A set 
of internal flat fields and arc frames were obtained immediately after each set of observations of 
W40 for flat fielding and wavelength calibration purposes. 
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Fig. 1. — 2MASS K-band image of the W40 IR stellar cluster, approximately 6' field of view (North 
is up, East to the left). Source designations are from Smith et al. (1985) except for IRS 5 which has 
not been previously reported. Many of the original IR sources reported by Smith et al. (1985) are 
revealed to be multiples. (This mosaic was created using Montage v3.0: Montage is maintained by 
the NASA Infrared Science Archive (IRS A), part of the Infrared Processing and Analysis Center 
(IPAC)). 



Table 1. Sources Observed with SpeX 



Source 


2MASS 


2MASS 


2MASS 








Exp. Time (s) 


Exp. Time (s) 




ID 


RA (J2000) 


Dec (J2000) 


J 


II 


K 


(SXD) 


(LXD) 


IRS 1A N a 


18312782-0205228 


18h31m27.81s 


-02d05m21.87s 


9.31 a 


7.49 a 


6.13 a 


900 b 


300 c 


IRS IAS" 


18312782-0205228 


18h31m27.83s 


-02d05m23.61s 


8.11" 


7.00 a 


6.49 a 


480 


300 c 


IRS IB 


18312866-0205297 


18h31m28.66s 


-02d05m29.78s 


11.47 


9.15 


7.50 


480 


120 


IRS 1C 


18312601-0205169 


18h31m26.02s 


-02d05ml6.97s 


10.87 


9.11 


8.08 


480 




IRS 2A 


18312397-0205295 


18h31m23.97s 


-02d05m29.51s 


8.74 


7.19 


6.00 


24(1 


120 


IRS 2B 


18312257-0205315 


18h31m22.58s 


-02d05m31.60s 


9.72 


8.79 


8.29 


480 




IRS 3A 


18312395-0204107 


18h31m23.95s 


-02d04ml0.76s 


8.91 


8.02 


7.46 


240 


360 


IRS 5 


18311482-0203497 


18h31ml4.82s 


-02d03m49.80s 


8.27 


7.51 


7.08 


240 





Note. — All SpcX observations carried out 08 July 2006 (UT) except where noted. JHK photometry is from 2MASS except where 
noted; relative uncertainty for the 2MASS JHK photometry is less than 2 %. 

a IRS 1A is not resolved in the 2MASS images. Positions listed for IRS 1A North and South are derived from our SpcX Guide Dog 
images; photometry is derived from the SpcX JHK spectra with a relative uncertainty of roughly 10 %. 

b Observed 09 July 2010 (UT). 

c Observed 13 May 2010 (UT). 



* 




Corrpanion? 

Fig. 2. — (a) SpeX guide camera K-band image of the region surrounding IRS 1 and 2. (b) Detail of 
IRS A showing resolution of 1A into North and South components, (c ) Detail of IRS 2A showing 
possible companion to the east. 
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In order to improve our S/N, we re-observed IRS 1A North and South on 13 May and 08-09 
July 2010. Conditions for the May observations were good: clear skies and ~ Qf!7 seeing in K. In 
July conditions were excellent, with photometric skies and seeing from (f.'3 - Of! 7 in K. Observational 
procedures were the same as for the 2006 observations except with longer integration times and 
additional pair cycles to improve S/N. Calibration and reductions were also identical to the 2006 
observation run. Note, however, that the S/N for IRS 1A North in the LXD mode was limited by 
the calibration source which was significantly fainter than IRS 1A North. 

All data were reduced using Spextool ( Cushing et al.pOOl ), the IDL-based package developed 
for the reduction of SpeX data. The Spextool package performs non-linearity corrections, flat field- 
ing, image pair subtraction, aperture definition, optimal extraction, and wavelength calibration. 
The sets of spectra resulting from the individual exposures were median combined and then cor- 
rected for telluric absorption and flux calibrated using the extracted AO V telluric standard spectra 
and the technique and software described by Vacca et al. ( 2003 ) . Uncertainty in the SpeX flux cali- 
bration is ~ 10 % ( Cushing et al.||2005 Rayner et al.||2009 ). The spectra from the individual orders 
were then spliced together by matching the flux levels in the overlapping wavelength regions, and 
regions of poor atmospheric transmission were removed. The complete, flux calibrated 0.8 — 5.0 /im 
spectra are shown in Figures [3] and |4j 



3.2. MIRSI Observations 

Mid-IR observations of the central portion of the W40 IR cluster were carried out at the IRTF 
using the Mid-InfraRed Spectrometer and Imager (MIRSI, Kassis et al.||2008 ) instrument on 8 July 



2006 (UT). MIRSI offers complete spectral coverage over the atmospheric windows at 8-14 /im and 
18-26 for both imaging (discrete filters and a circular variable filter) and spectroscopy (in the 
10 and 20 //m windows with resolutions of 200 and 100, respectively) utilizing a Raytheon 320x240 
pixel Si: As array with a plate scale of 0'. / 27 per pixel, and a field-of-view of 86"x63" at the IRTF. 

Images of the central 3' of the IR cluster were obtained using the 8.7 and 12.3 p filters. 
Approximately 15 pointings were used to cover IRS 1 through 3 and IRS 5 with on-source exposure 
times of 25 s for each pointing in each filter. Standard chop and nod techniques were used for 
each pointing to eliminate infrared background emission. Standard stars a Lyr (A0V) and fi Cep 
(M2Ia) were observed throughout the night for flux calibration. The final mosaics are shown in 
Figure [5] Note that only the region surrounding IRS 1 - 3 is shown as IRS 5 was not detected 
in any of the filters used. The infrared sources IE and F, and 2 D - F are reported here for the 



first time and have been named according to the convention of Smith et al. (1985). Fluxes derived 
using standard aperture photometry techniques are shown in Table [2} For the brighter sources, 
the uncertainty is dominated by systematic errors in the flux calibration, whereas for the fainter 
sources, it is dominated by uncertainty in the background emission. 
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Fig. 3. — Spectral energy distributions for main sequence sources in W40. Solid line — SpeX near-IR 
spectra (this study); Crosses — MIRSI mid-IR photometry (this study); Triangles — Spitzer-IRAC 
photometry (this study); Diamonds — Smith et al. (1985) photometry; Dotted line — stellar atmo- 



spheric model fit for given reddening and distance (see text). 




Fig. 4. — Spectral energy distributions for young stellar objects and pre-main sequence stars in W40. 
Line — SpeX near-IR spectra (this study); Crosses — MIRSI photometry (this study); Triangles — 
Spitzer-IRAC photometry (this study); Diamonds — Smith et al. (1985) photometry. 
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Fig. 5.— MIRSI mosaic of W40 IRS 1 - 3 at 8.7 and 12.3 /an. Final mosaic is 2.6' x 2.3 / (North is 
up, East is left). 



Table 2. MIRSI photometry 



Source 


8.7 /zm 


12.3 fim 


Notes 




mjy (rel. err) 


mjy (rel. err) 




IRS 1A 


5040.6 (0.06) 


4104.2 (0.06) 




IRS IB 


2115.0 (0.1) 


2602.7 (0.1) 




IRS 1C 


381.3 (0.1) 


468.5 (0.1) 




IRS ID 


176.6 (0.1) 


<400 




IRS IE 


289.8 (0.1) 


240.8 (0.1) 




IRS IF 


364.9 (0.1) 


1008.8 (0.1) 


Diffuse 


IRS 2A 


28571.5 (0.06) 


43307.4 (0.06) 


17140.2 (0.2) 


IRS 2B 


357.1 (0.1) 


<400 




IRS 2C 


<300 


<400 


Not detected 


IRS 2D 


99.1 (0.1) 


<400 




IRS 2E 


255.7 (0.1) 


<400 




IRS 2F 


394.4 (0.1) 


241.1 (0.1) 


Diffuse 


IRS 3A 


741.9 (0.1) 


1952.4 (0.1) 


Resolved, non-point-likc 


IRS 5 


<300 


<400 


Not detected 



Note. — An aperture radius of 3" was used for all measurements. 
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4. Spitzer Space Telescope Archive Data 

We also retrieved IRAC images of the W40 central cluster from the Spitzer Space Telescope 
Archive, which were originally obtained as part of Spitzer Program P30576 (Allen et al.). These 
observations were carried out in IRAC "High Dynamic Range" (HDR) mode so that a short (0.6 s) 
and long (12.0 s) exposure was obtained for each discrete pointing. The central part of W40 suffers 
from numerous saturation effects/artifacts in the long exposure data, so we restricted our study to 
only the short exposure images. We obtained the "basic calibrated data" (BCD) for all IRAC bands 
within the central 10' of IRS 2A using the Spitzer Heritage Archiv^] and then assembled a mosaic 
for each band (1, 2, 3, and 4) using MOPEX0 ( |Makovoz & Khan|[2005| with default settings. A 



three-color composite image at 3.6, 5.8, and 8.0 jim. (channels 1, 3, and 4) is shown in Figure |6j 
Simple aperture photometry was carried out to determine the fluxes of the brightest sources in 
each channel using the prescriptions and aperture corrections in the IRAC Instrument Handbook^ 
For all targets we picked the largest aperture possible that also did not contain any contaminating 
flux from the PRF wings of nearby sources. The appropriate aperture corrections were applied to 
all measurements. For IRS5 in Band 4, the background nebular emission varies substantially over 
small scales close the the star: we chose a smaller aperture in this case to minimize contamination 
from the background nebula. No color correction has been applied. Fluxes and aperture sizes are 
listed in Table [3j According to the Handbook, errors should be <10% for channels 2 — 4 and < 
15% for channel 1. A few of the brightest sources are still saturated in some of the bands. Since 



IRS 3 A does not appear to be a point source (see Sect. 6.2), a large (12") aperture was used with 
no correction. 



5. Analysis & Results 

5.1. Spectral Energy Distributions 

In Figures [3] and [4] we plot the full spectral energy distributions for each source from the optical 
through the mid-IR, including data obtained in this study (SpeX, MIRSI, and Spitzer Archive) and 



data from Smith et al. (1985). The agreement among the datasets at overlapping wavelengths is 



generally quite good with two notable exceptions: 



• IRS1A: The Smith et al. (1985), MIRSI, and Spitzer images do not resolve IRS 1A into its 



northern and southern components, as seen in our SpeX images (see Section 3.1 ). Comparing 



the near-IR spectra for N and S to the combined fluxes reported in Smith et al. (1985) and 



1 


http : // sha. : 


ipac . caltech . edu/applications/Spitzer/SHA/ 




2 


http : // irsa 


ipac . caltech. edu/data/SPITZER/docs/dataanalysistools/tools/mopex/ 


3 


http : // irsa 


ipac . caltech. edu/data/SPITZER/docs/irac/iracinstrumenthandbook/ 
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Fig. 6. — Spitzer-IRAC composite mosaic of W40 at 8.0 (red), 5.8 (green), and 3.6 jum (blue). 
Field-of-view is roughly 10' (North is up, East is left). 
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our MIRSI data, it appears that the longer wavelength emission (>4 //m) is dominated by 
IRS 1A North, whereas the shorter wavelengths are dominated by IRS 1A South (Fig. [7]). In 
the separate SEDs for the northern and southern components, we divide up the fluxes from 



Smith et al. (1985) and MIRSI by assigning all data shorter than 1.5 /im to IRS 1A South, 



and all points longer than 4 /im to IRS 1A North. 

IRS3A: The absolute levels of the SpeX spectra in the J, H, and K bands do not agree well 



with either the long wavelength (LXD) SpeX data or the near-IR photometry from Smith 



et al. (1985). Since the long and short wavelength SpeX data were obtained on the same 
night, it seems unlikely that this is due to intrinsic variability and instead is due to some 
issue with the observations. Indeed the first part of the night during which 3A was observed 
suffered from variable seeing. In this case, we decided to scale the J, H, and K band SpeX 



data by a factor of 0.75 to match the Smith et al. (1985) fluxes (scaled fluxes for IRS 3A are 



shown in Fig. |3j). IRS 3A also has notable excess flux in the 8.0 /im IRAC band (channel 
4) and the 12.3 /jm MIRSI band. The excess at 8.0 /im is likely due to enhanced emission 
from the PAH band at 7.7 /im. The excess at 12.3 /im, however, could be due to the 12.7 /im 
PAH line, crystalline silicates, 12.82 /jm [Ne II] emission, or some combination of all three. 
(Since the bandpass for the 12.3 /im filter is only 9.6%, emission from the strong PAH feature 
at 11.3 /im would not be included.) Further mid-IR spectroscopy of IRS 3A is required to 
discern the source of the excesses. 

Except for IRS 5, all the sources show significant emission in the 5-12 /im range indicative 
of heated circumstellar dust. The mid-IR spectral index, a m i r = dlog(\F\)/d\, for each source 
was determined by fitting a line to the \F\ values between 3 and 10 /im. The results are listed in 
Table H 



5.2. Near-IR Spectra: Main Sequence Stars 



Four of the observed sources — IRS 1A South, 2B, 3A, and 5 — have strong Brackett and Paschen 
series absorption indicative of late-0 and early-B type stars (Fig. [8]). None of these sources has 
obvious C IV emission (2.078 /im), N III emission (2.1155 /im), or He II absorption (2.1885 /im) 
which restricts the spectral types to 09 or later. Determining spectral types based on line equivalent 
widths alone is notoriously difficult ( Hanson et al.|[2~005| ), so we have compared our spectra directly 
with early-type stars from the IRTF spectral library (Vacca et. al., in preparation). Derived spectral 
types are listed in Table [4] and are accurate to about +/- 1 sub- type. We also measured the 
equivalent widths of He I (2.1126 /im) and Br7 (2.16 /im) lines for each source and compared them 
to the Hanson et al. (1996]) and |Blum et al. (1997) spectral atlases. In all cases, our EWs are 
consistent with early-B spectral types, but note that it is very difficult to discriminate between 
subtypes based on EWs alone. The existence of O and B main sequence stars at the center of the 
W40 IR cluster implies an upper limit on the age of the W40 cluster of ~ 7 My. 



-14- 



Table 3. IRAC Photometry 



Source 


3.6 (im a 


4.5 /W 3 


5.8 nm b 


8.0 [im b 


Aperture 




mJy 


mJy 


mJy 


mjy 


Radius (") c 


IRS 1A 


>950 


>980 


>6950 


>3700 


SAT 


IRS IB 


>950 


>980 


1742.4 


>3700 


6.0 


IRS 1C 


445.1 


410.1 


348.5 


301.0 


3.6 


IRS 2A 


>950 


>980 


>6950 


>3700 


SAT 


IRS 2B 


192.0 


146.9 


153.6 


263.0 


3.6 


IRS 3A 


452.1 


452.0 


535.5 


1765.0 


12 c 


IRS 5 


535.8 


410.6 


410.0 


311.8 


6 (3.6) d 



a Relative error: 15% 
b Relative error: 10% 

c DifTerent aperture radii were used to avoid contamination from 
nearby soruces; see text for discussion. 

d An aperture radius of 3'.'6 was used for the 8.0 fim image to avoid 
contamination from variations in the nebular background close to the 
star. 



° Since IRS 3 A does not appear to be a point source (see Sect. [672 
a large aperture was used with no correction. 



RS 1A S (09.5) , 



E 10 - ' 




RS 1A N (HAeB 



Log A (jj.rr 



Fig. 7. — Spectral energy distribution for IRS la. Solid line — SpeX near-IR spectra (this study); 
Grey line — combined near-IR spectrum of 1A North and South; Crosses — MIRSI mid-IR photom- 
etry (this study); Diamonds — Smith et al. (1985) photometry. See text for discussion. 
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09.5 



1 .5 
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•Si 

b 

+ 



X 
=5 



IRS 5 

B1 



E 1 .0 

o 



RS 3A 
B3 



IRS 2B 



0.5 



B4 



J I I L 



_L 



J I I L 



J_ 



y^^yA 




J I I L 



J_ 



J I I l_ 



1 .0 



.5 



2.0 



Wavelength (/U.m) 

Fig. 8. — Normalized near-IR spectra for candidate Main Sequence stars in our sample. 
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All of the main sequence sources in our sample exhibit deep, unresolved He I absorption at 
1.083 fxra (see Table [5]), which is often the deepest line in the spectrum. Though relatively rare for 
"normal" O and B-type stars ( Conti fc Howarth||1999 ), it is quite prominent in both emission and 
absorption for massive stars with dense winds including some Be stars QGroh et al.|2007| ). Since the 
line widths are very narrow, and none of our main sequence (MS) stars exhibit Be star properties, 
it seems unlikely that the strong He I absorption is due to stellar winds. Note also that none of 
our MS stars, except IRS 5, is a strong X-ray emitter, which is also a signature of stellar winds in 



B stars (Cassinelli et al. 1994). 



Relatively strong He I 1.083 /im absorption can also be produced by non-LTE effects in the 



atmospheres of early- type stars (Lennon & Dufton 1989; Przybilla 2005). However, the equivalent 



widths we measure for the W40 sources are far larger than those measured by Lennon & Dufton 



(1989), or derived from the models of Przybilla (2005). In fact, the 1.083 /j,m absorption is also 



far stronger than that observed in the reference spectra of the IRTF spectral library that were 
used to determine the spectral types of the W40 sources. These results strongly suggest that the 
1.083 /im absorption is not due to non-LTE effects in the stellar atmospheres of these sources. He I 
absorption is also thought to be associated with coronal magnetic activity, and hence is commonly 



observed for low- mass main sequence stars and close binaries (Zirin 1982). It seems likely that all 



of our MS sources are in fact binaries with close and/or low-mass companions. In addition, IRS 5 
is a known X-ray source ( Kuhn et al.||2010 ) which further supports this conclusion for that source. 



For each main sequence star in our small sample, we simultaneously determined the extinction 
(Ay) and distance by fitting a reddened atmospheric model to the SpeX near-IR spectra using a x 2 



minimization routine. We adopted a Kurucz standard model atmosphere (Kurucz 1993) for each 



target with the stellar radius and effective temperature for each spectral type taken from|de Jager &: 



Nieuwenhuijzen (1987). We adopted the reddening law from Fitzpatrick & Massa (2007) for a fixed 



Ry of 3.1. The results, including our adopted radii and and effective temperatures, are presented 
in Table [5] with resulting reddened model atmospheres plotted in Figure [3j In general, the fits to 
the near-IR data are excellent; the values for distance and Ay agree with previous studies quite 
well. Note that the reddened model SED for both IRS 1A South and IRS 3A is always well below 



the UBV photometry from Smith et al. (1985), possibly due to uncertainty in the optical-to-IR 
extinction law for W40. 

The dominant source of error in the fits is the ±1 subtype uncertainty in the spectral type 
whereas the dominant source of error in Ay is the extinction law. To estimate the uncertainties 
(particularly for the distance), we ran the fitting routine for +/- 1 spectral-subtype (with an 
associated change in R and T e jf) for each source with the resulting range in distance shown in 
Table [5j For IRS 1A South, we also tried R and T e fj values from Vacca et al. (1996) and Martins 
et al.| ( |2005 ). The Vacca et al. (1996) values for R and T e fj are somewhat larger than those of 
de Jager & Nieuwenhuijzen (1987), which pushes IRS 1A South about 100 pc further away. The 



value for T e ff from Martins et al. (2005) is roughly the same as that of de Jager & Nieuwenhuijzen 



(1987), but Martins et al. give a much smaller radius for an 09.5 star, which brings IRS 1A South 
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Table 4. Spectral Features and Source Properties 



Source 


Sp. Type 


He I (1.083 /Ltm) 
EW (A) 




CXOW40? a 


3.6 cm? b 


IRS 1A North 


HAe/Bc 


abs/ cm 


-1.3±0.2 


Yes? 


Yes 


IRS 1A South 


09.5 


1.81 




Yes? 


No 


IRS IB 


Class II 


P Cygni 


-0.6 ±0.1 


Yes 


Yes 


IRS 1C 


Class II 


em 


-0.9 ±0.1 


Yes 


Yes 


IRS 2A 


HAe/Be 


P Cygni 


1.0 ±0.1 


Yes 


Yes 


IRS 2B 


B4 


1.92 


-0.2 ±0.3 


No 


No 


IRS 3A 


B3 


1.90 


-0.3 ±0.1 


No 


N/A 


IRS 5 


Bl 


1.25 


-1.6 ±0.1 


Yes 


N/A 



^Kuhn et aT] pOlO] 



Rodriguez et al. (2010 i 



Table 5. Summary of Near-IR Atmospheric Model Fits to Main Sequence Sources 



Source 


Sp. Type 




Radius a 


Model T el f 


log(Qo) 


E(B - V) h 


A v 


D 


D - range 






(K) 


(«©) 


(K) 








(pc) 


(PC) 


IRS 1A South 


09.5 


32210 


7.71 


32000 


47.8 


3.4 


10.6 


536.0 


441—578 


IRS 2B 


B4 


16904 


3.43 


17000 


42.8 


2.9 


9.1 


455.5 


396—526 


IRS 3A 


B3 


18793 


3.76 


19000 


43.7 


3.0 


9.2 


349.0 


312—416 


IRS 5 


Bl 


25410 


5.16 


25000 


45.7 


2.8 


8.5 


469.1 


340—686 



^de Jager &: Nie"uw cnhuijzcn| ( |1987| > 



b All fits assume Ry = 3.1. 
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about 40 pc closer. 



Although the acceptable ranges in distance for each source overlap very nicely, there are some 
discrepancies among the best-fit values derived for each source. The agreement in distance between 
IRS 2B and 5 is striking (< 3%, well within the estimated errors). The large error in distance for 
IRS 5 is due to the large spread in effective temperature for early B-stars. IRS 3A appears to be 
much closer than the others which does not seem physically plausible given the nearly identical Ay 
and assumption of cluster membership. In addition, the model fit for this source is not as good as 
the fits for the other sources. This may be due to a possible close binary companion, as suggested 
above (IRS 3A will be discussed further in Section 6.2). On the other hand, IRS 1A South appears 
to be somewhat farther than both IRS 2B and 5. This could be rectified if 1A South were one 
subtype later and IRS 2B and 5 roughly one subtype earlier. An alternative explanation is that 
the distance to 1A South is correct, and that both IRS 2B and 5 are overluminous. This would 
further support the argument that these two objects are in fact close binaries, as indicated by the 
strong He I 1.083 /im absorption. Changes in the assumed extinction law do not affect the distance 
strongly, but can produce different values for Ay. Changing the spectral type by one sub- type (as 
above) generally causes a variation in Ay of less than a few percent. Adopting a value of Ry of 
5.0 causes a decrease in the Ay of roughly 20%. 

Except for IRS 1A South (which is not resolved from 1A North in the mid-IR), all of the 
main sequence objects show at least some mid-IR excess indicative of heated circumstellar dust. 
Based on the combined fluxes (Section 5.1), we have attributed the majority of the mid-IR flux 
from IRS 1A to the northern component. It is, however, possible that some of the mid-IR flux 
is due to circumstellar dust associated with IRS 1A South. For IRS 2B, the IR excess does not 
become substantial until ~ 6 ^m, and is very strong at 10 /um. It is not clear whether the flux 
at 10 /iin is due to silicate or thermal emission; further observations at longer wavelengths are 
required. In either case, a flux deficit in the near-IR with strong IR emission at longer wavelengths 



is often modelled by an optically thick disk with a central "hole" (Calvet et al. 112005 D'Alessio et al. 



2005), which may be created by photoevaporation (Clarke et al. 2001) or the presence of a close 



companion ( Uchida et al.|2004 Rice et al.|2003 ), or both. On the other hand, the mid-IR excess for 
IRS 5 decreases uniformly, suggesting a much cooler, older disk. The mid-IR excess for IRS 3A is 
very flat from 3 to 10 ^m punctuated by excess emission in the 8 and 12 /im bands; this source will 
be discussed further in Section 6.2 The presence of circumstellar dust for these sources indicates 
that there has been insufficient time for these stars to dissipate their natal disks, and hence they 
may have only recently arrived on the main sequence. 



5.3. Near-IR Spectra: Young Stellar Objects 

The sources IRS 1A North, IRS IB, IRS 1C, and IRS 2A all have very different near-IR spectral 
features compared to the Main Sequence stars identified in our sample (see Fig. [9] and Table [4|. 
Both IRS 1A North and IRS 2A have strong Brackett and Paschen series emission, coupled with 
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strong absorption in the higher level Paschen lines. Both stars show Ca II and O I in emission 
as well. Based on these spectral features, and the strong mid-IR emission for both, we conclude 
that IRS 1A North and IRS 2A are likely heavily reddened Herbig Be stars. Using the strongest 
hydrogen emission lines in both sources, we derive an Ay of « 10, assuming Case B recombination, 
which agrees with that derived from the main sequence continuum fitting. All of the YSOs in our 
sample exhibit strong 1.083 [im. emission and absorption (including P Cygni profiles for IRS IB 



and IRS 2 A) which is common in both T Tauri and HAeBe stars (Edwards et al. 1987). 



The near-IR spectra for IRS IB and 1C are generally smooth except for modest Brackett 
and Paschen series emission coupled with strong Na I, Ca I, and CO absorption. The near-IR 
spectral features and the mid-IR spectral index (a m i r ) suggest that both of these sources are 
heavily reddened Class II low mass Young Stellar Objects (YSOs) ( [Greene & Ladajl996l >. 
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Fig. 9. — Final SpeX flux-calibrated spectra: Pre-Main Sequence and Young Stellar Objects. 
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6. Discussion 



6.1. The Distance to W40 



Depending on the technique used, the distance to the W40 complex has been estimated to 



be between 300 and 900 pc (Rodney & Reipurth 2008; Kuhn et al. 2010 and references therein). 
If we exclude IRS 3A (which has a relatively poor fit, and is also exceptional in other ways, see 



Section 6.2), then the best estimate for the distance based on our results (for 3 stars) is between 
455 pc and 536 pc. We can confidently rule out distances less than 340 pc and greater than 686 pc; 
making it very unlikely that W40 is part of the same complex that includes the Serpens star forming 



regions, as assumed by Bontemps et al. (2010). Our best estimates are somewhat lower than the 



distance estimate of 600 pc derived from X-ray luminosity function (XLF) fitting method used 



by Kuhn et al. (2010). (In fact, a distance of 500 pc produces a qualitatively worse fit to the 
XLF according to the authors.) Our distance determination has the advantage of being based on 
(admittedly only a few) well-determined spectral types, but only for a few sources that happen to 
be OB stars with a wide range of possible luminosities. The advantage of the XLF method is that 
it is based on statistics over a large number of observed sources; the disadvantage is that one must 
assume a "universal XLF" that may not apply to all star forming regions. For the remainder of 
the paper, we will adopt a distance of 500 pc. 



6.2. IRS 3A 



As indicated in Section 5.1 IRS 3 A has significant mid-IR excess with a generally flat spectral 
index (a m i r = —0.3). In addition, IRS 3A is also resolved as a circular source at longer mid-IR 
wavelengths. This can be seen readily in both the MIRSI 8.7 and 12.3 /im images (Fig. [5), and m 
the IRAC channel 4 images (see red channel in Fig. [6]). In fact, IRS 3 A is point-like only in the 
IRAC channel 1 (3.6 /im) data; in the other bands it has a roughly Gaussian profile with a FWHM 
that exceeds that of the corresponding PSF. The measured FWHM as a function of wavelength 
(with comparison to the PSF FWHM) is shown in Table [6} At 500 pc, the measured FWHM at 
4.5 /im corresponds to a size of 1050 AU, increasing to 2700 AU at 12.3 /jm. The fact that IRS 3A 
is resolved in the near- to mid-IR indicates that the emission must be coming from a circularly 
symmetric distribution of dust — either a spherical envelope or a disk viewed roughly face-on — 
with a radius of roughly 500 AU. There are, however, issues with both of these interpretations. 
In general, near- to mid-IR emission from circumstellar disks around intermediate mass stars (e.g. 
HAeBe stars) is generated at the disk surface where dust temperatures exceed 250 K. For most 
flaring disk models, this corresponds to regions at the disk surface within 10 AU of the central 
star (e.g. Dullemond||2002 ) , and hence would be unresolved in our observations. On the other hand, 



an optically thick spherical envelope would produce additional extinction in the optical that would 
produce a larger Ay for IRS 3A than the other stars in W40, which is not observed. A more 
probable explanation is that the resolved emission from 4-12 fim is due to emission from small 
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grains (transiently heated by the strong UV emission from IRS 3A) in an optically thin, dusty 



nebula (or shell) immediately surrounding the star (e.g. Hartmann et al. 1993). If this is the case, 
then PAH and [Ne II] emission from this dusty nebula could also explain the excess flux observed in 
the 8.0 and 12.3 fj,m bands (see Figure [3]). This dusty shell may be a relic from an earlier pre-main 
sequence (Herbig Be) phase which has not yet been dispersed by the stellar winds of IRS 3 A. 

The fact that the Kurucz model fits yield a distance to W40 that is far closer than that of 
the other MS stars in our (admittedly small) sample, and the presence of strong He I absorption 
strongly suggests that IRS 3A is in fact a close binary. Increasing the distance of IRS 3A by a 
factor ~1.3 would place the binary at the same distance as IRS 2B and IRS 5, and would result 
in a ~1.7 fold increase in the luminosity, which could be accounted for by a companion of mid-B 
spectral type. It is also possible that IRS 3A is overluminous because it is a "bloated" star, i.e. the 
stellar radius is much larger than a typical MS B3 (e.g Ochsendorf et aL||2011 ). Pre-main sequence 
evolutionary tracks for intermediate mass stars show a significant increase in total luminosity just 



before the star settles onto the main sequence (Palla Sz Stahler 1993; Hosokawa & Omukai 2009). 



It is possible that IRS 3A is nearing (or at) this peak in luminosity, which would contribute to 
the poor stellar model fit, result in a shortened distance, and explain the existence of the extended 
dusty envelope. In this case, the strong He I absorption could still be due to a much lower mass 
binary companion. 



6.3. Compact Radio Emission 



All of the sources we identify as YSOs in our sample are also detected in the radio survey 
of Rodriguez et al. (2010) (see Table [4]). Only two of the main sequence stars in our study are also 



included by the radio maps of Rodriguez et al. (2010) (IRS 1A South, and IRS 2B) and neither is 



detected at 3.6 cm. Rodriguez et al. (2010) suggest that the non- variable compact radio sources 



are likely due to ultra-compact H II (UCH II) regions around young, high mass stars. This seems 
unlikely, however, as the size of the unresolved radio sources at a distance of 500 pc would be less 
than ~ 100 AU — much smaller than a typical UCH II regionf Kurtz 2005). It seems more likely 



that the 3.6 cm continuum flux would be due to free- free emission from shocked gas within 100 AU 
of the YSO caused by a jet or outflow. It could also be driven by synchrotron emission associated 
with the stellar magnetic field, a conclusion that is supported by the fact that all the YSOs in our 
sample are also strong X-ray sources (see Table El) . 



6.4. The W40 H II Region 



The predicted Lyman continuum (LyC) fluxes for each of our main sequence sources (based on 
the Kurucz (1993) model) is given in Table [5| The total LyC flux in the region is dominated by the 
09.5 star IRS 1A South at 6.3xl0 47 photons s _1 . Based on our spectral type analysis (Sect |5.2[), 
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it is unlikely that IRS 1A South is earlier than 09 which places an upper limit on the LyC flux 
available to drive the H II region of ~ 1.6 x 10 48 photons s . This upper limit agrees very nicely 
with the derived LyC flux of ~ 1.5 x 10 48 photons s~ x from radio observations (Goss & Shaver 



1970 Smith et al.||1985 ). Based on the IR luminosity, Smith et al. (1985) originally suggested that 
IRS 2A must be the dominant energy source in W40. We conclude, however, that IRS 1A South 
must be the main source of ionizing radiation. 

The morphoplogy of the H II region in the thermal IR is quite different from that seen in the 
POSS plates and 2MASS images. Images of the region from the Mid-course Space Experiment 
(MSX) obtained from the NASA Infrared Science Archive (IRSA)Q (see Fig. [To]) reveal an hour- 
glass shaped structure centered roughly on IRS 2 A, approximately 28' in height and about 17' in 
width (Rodney & Reipurth 2008), corresponding to 4.1 by 2.5 pc at a distance of 500 pc. The 
long axis of the hourglass is at a PA of ~ 300 degrees (NNW). The structure is most prominent 
in the MSX A-band (8.3 [im.) images, where the brightness is likely dominated by emission in 
the PAH line at 7.7 /im. The hourglass does not appear related to the weak CO outflow centered 



on the dense material roughly 3' to the West (Zhu et al. 2006). We interpret this double- lobed 



morphology as a pinched waist "bubble" blown by a symmetric wind, most likely from the 09.5 star 
IRS 1A South. The waist axis of the hourglass lies along a PA of roughly 60 degrees (ENE). The 
Spitzer IRAC images (Fig. [6]) reveal numerous bright-rimmed sheets and clumps ("elephant trunks") 
along this axis, nearly all of which point back toward IRS 1A, further underscoring IRS 1A South 
as the dominant wind and UV source of the region. Far-IR and mm-wave observations of W40 
reveal 36 deeply embedded sources (prestellar cores, and Class 0,1 objects) indicating that star 
formation is on-going, and has perhaps been initiated by the pressure wave associated with the 
wind-bubble ( |Maury et al]|2011[ ). 



http : // irsa. ipac . caltech. edu/Missions/msx .html 
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Table 6. IRS 3A: Size vs. Wavelength 



Waveband 


PSF a 


IRS 3A 


IRS 3A 


(H 


FWHM (") 


FWHM (") 


FWHM (AU) b 


IRACl, 3.6 


1.66 


1.764 


882 


IRAC2, 4.5 


1.72 


2.1 


1050 


IRAC3, 5.8 


1.88 


2.58 


1290 


IRAC4, 8.0 


1.98 


4.158 


2079 


MIRSI, 8.7 


1.21 


4.32 


2160 


MIRSI, 12.3 


1.32 


5.4 


2700 



a IRAC PSF taken from mean FWHM of the point response 
function (PRF) reported in the IRAC Instrument Handbook. 
MIRSI PSF based on measured performance and is nearly 
diffraction limited. 

b Assuming d = 500 pc. 




Fig. 10.— MSX 3-color composite of W40 at 8.3 (blue), 12.1 (green), and 14.6 fim (red). Field of 
view is approximately 30'. Red circle marks the center of the peak CO emission (Zhu et al. 2006). 
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7. Summary 

We have presented the first detailed analysis of the eight brightest members of the galactic H II 
region W40 based on medium resolution near-IR spectroscopy (IRTF/SpeX) and mid-IR imaging 
(IRTF/MIRSI, Spitzer IRAC). SpeX guide camera images reveal that IRS 1A is in fact a binary 
(with "North" and "South" components). In addition, IRS ID breaks up into at least 7 distinct 
sources in the near-IR. Our MIRSI images reveal 5 new sources in the mid-IR (IRS 1 E and F, 
IRS 2 D - F), two of which are diffuse in nature. 

Using the SpeX near-IR spectra we have classified all 8 of the sources observed based on near- 
IR emission and absorption features (see Table [4]). Four of the objects appear to be late-0 or 
early-B stars. Two of the stars have spectral features indicative of Herbig Ae/Be stars, and the 
remaining two appear to be low mass YSOs (Class II). Strong He I absorption at 1.083 /xm strongly 
suggests that the main sequence objects may in fact be close binaries. Spectral energy distributions 
reveal that 7 of the 8 sources have mid-IR excesses indicative of circumstellar disks. 

For the main sequence sources, we simultaneously derived Ay and D for each star by fitting a 
reddened stellar atmosphere appropriate to the spectral type. In general, the agreement between 
the sources is good. We find a best estimate for the distance of 455 — 535 pc with Ay ranging 
from 8.5 to 10.6. Our derived distance agrees well with previous studies, though is somewhat 



lower than the recent determination by Kuhn et al. (2010) based on XLF fitting. Our distance 



determination indicates that the radio emission from some of the sources reported in Rodriguez 



et al. (2010) cannot be due to UCH II regions (due to size constraintss) and is instead most likely 



due to shocked gas due to jets and/or outflows from the source. 

IRS 3A is found to be somewhat peculiar in a few ways: first, despite having a fairly clear B3 
spectral type, stellar model fits were quite poor and the distance to the source was found to be 
significantly smaller than that for the 3 other MS stars. Furthermore, IRS 3A appears to have a 
substantial dusty envelope, at least 2700 AU in size, that is resolved at mid-IR wavelengths. We 
suggest that IRS 3A is either a close binary composed of two early-B stars, or a PMS star with a 
close companion and/or a strong stellar wind. 

Our observations indicate that the 09.5 star IRS 1A South is the dominatant source of LyC 



photons for the W40 H II region (as opposed to IRS 2A, as originally suggested by Smith et al. 



(1985)). Furthermore, IRS 1A South is situated at the center of the pinched waist morphology 
seen in wide-field mid-IR images (MSX), suggesting that it is the source of the stellar wind that 
has blown this bubble over time. 



This research has been supported in part by the Universities Space Research Assoc. (USRA) 
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